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At the end of their lives massive stars can form neutron stars in their womb, when they 
explode as core-collapse supemovae. Neutron stars were proposed theoretically by Baade 
and Zwicky in 1934, but the observational confirmation of their existence only came in 
1967 when they were discovered by Bell and Hewish as pulsars. The interpretation of 
these pulsars as rapidly rotating magnetized neutron stars was established by 0]|. This 
type of pulsars, in which the rotation of the neutron star is responsible for the observed 
luminosity, are known as Rotation-Powered Pulsars (RPP). The spin evolution of RPP (in 
the canonical model) is well known and used to estimate ages and surface magnetic fields of 
neutron stars (see Fig. |TJ). Here, we briefly summarize the main properties of the canonical 
model, in order to do a comparative analysis with our results. 

The energy loss by radiation allows the neutron star to undergo a systematic spin-down 
(El). If a pulsar spins down from an initial spin period F\ h then, the deceleration of the 
pulsar, Q, is given by an empirical formula obtained by balancing the spin-down luminosity 
with the energy loss by radiation (a dipole magnetic field), 
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2.51(1) 

-1.5(1) 

2.140(9) 

1.4(2) 

2.684(2) 

2.839(3) 

2.65(1) 

0.9(2) 

Lyne 1988 
Middleditch 2006 
Nagase 1990 

Lyne 1996 
Weltevrede 2011 
Kaspi 1994 
Livingstone 2007 
Espinoza 2011 


Table 1: Observational datas of young pulsars (t < 10 kyr). The periods and period deriva¬ 
tives are taken from Manchester 2005. The braking indices obtained from observations, 
n 0 b s , are shown with their respective uncertainty in the last digit. 


where k is a constant that accounts for structural factors of the pulsar, n is the braking 
index and Q = 2tt/P is angular velocity of the pulsar. The case n — 3 correspond to 
purely dipolar radiation. Different values of n would correspond to different processes of 
rotational energy loss. The characteristic age of the pulsar is given by, 
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Figure 1: Dipolar oblique rotator schematic: The pulsar rotate with angular velocity Q and 
it has a magnetic moment m inclinated an angle a respect to the rotation axis. 


n 

(n — l)f2 



P 

(n — l)P 



( 2 ) 


which is valid only when n^l. The parameter Q 0 is the initial angular velocity at t — 0. 
In the limit when f2 0 D (or equivalently P 0 -C P), the standard characteristic age of 
the pulsar is obtained: r = —Vt/(n — 1)D = P/(n — 1 )P = P/2P, assuming n = 3. 
Traditionally, this expression is taken as the definition of characteristic age, despite the fact 
that n can be different from 3, and is, in fact, so for all cases in which a stable, accurate 
value has been determined. In addition, if n is a constant of the pulsar then the spin-down 
luminosity and the spin period evolve with time according to, 
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where To = Po/{n — 1)P 0 is the initial spin-down timescale of the pulsar and E 0 is the 
initial spin-down luminosity, which has roughly constant energy output, E ~ E 0 , until a 
time r 0 beyond which evolve as E oc Similary, P ps P 0 for t -C r 0 and 

evolves at later times as P oc ([?]). Notice that changes in the moment of inertia 

are neglected as of now (we will discuss such possibilities in the second part of this paper). 
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The braking index n can be obtained directly from timing observations, and its measure¬ 
ments are crucial to understanding the physics behind the pulsar. Its definition is obtained 
from the observables as, n = QQ/Q 2 = 2 — PP/P 2 . The standard requirements for ac¬ 
curate measurements of n are: (a) that the pulsar is slowing fast enough to measure Q, 

(b) known position of the pulsar (~ 1”) and (c) that the braking is not affected by other 
agents such as glitches, timing noise or low frequencies. The glitches are important for 
t ~ 10 — 15 kyr. The timing noise varies from pulsar to pulsar, but it is correlated with 
the spin-down and can affect n if the observational data set is not very large ( [[3| ). The 
difficulty of extending this analysis to older pulsars lies in the fact that for these pulsars 
unfortunately measurements of Q require tens of years and the expected Q are very small. 
Furthermore for older pulsars, the aforementioned effects, such as glitches and timing noise 
are more recurrent which make such measurements more challenging (0|). Thus young 
pulsars (t < 10 kyr) are the best choice for measuring n. To date there have been only a 
few measurements for the braking index, and only for very few young pulsars. In all cases 
n 0 bs < 3 (see the Tableland references therein). For these young pulsar glitches are minor 
and the pulsar spin-down is faster. These results suggest that there should be present more 
complicated processes in the pulsar decreasing n from the expected value for the dipolar 
oblique rotator (n = 3). The standard factors that may affect the braking index are (CD): (a) 
multipolar electromagnetic radiation, n > 5, (b) quadrupole gravitational radiation, n — 5, 

(c) decay of the magnetic field, n > 3, (d) radial deformation of the magnetic field lines, 
1 < n < 3, (e) relativistic winds, n < 3, (f) transverse velocity of the pulsar, n < 3. 

More exotic factors are: (g) Intense emission of neutrinos in the early evolution of the 
pulsar, n < 0 (161), (h) Crustal movement of the neutron star by tectonic plates can produce 
n > 3 or n < 3 (0), (i) Growth of the magnetic field due to thermomagnetic instabilities 
in the crust of the neutron star, n < 3 ( Ii28l ). (]) Growth of an intense magnetic field 
submerged on neutron star crust in the hypercritical accretion phase, which re-emerge by 
ohmic diffusion, n < 3 m , m, ma, mox (k) Changes in the moment of inertia of the 
neutron star, n < 3 m, ca, my 

The effects of alignment ( da/dt > 0) or misalignment ( da/dt < 0) on n were calculated 
theoretically by [fI51 . For a dipolar field, n changes from its above canonical value of 3 
by an amount 2r (dinsin 2 a/di) due to alignment or misalignment between the rotational 
and magnetic axes, where r is the characteristic age introduced above. Nevertheless, ob¬ 
servations of alignment or misalignment are very complicated to carry out because large 
amounts of data are required to make a statistical estimate. In addition, unfortunately, an 
exact expression for the electromagnetic (Poynting) flux energy loss as a function of the 
magnetic inclination angle a remains still elusive. For these reasons, in this work we do 
not take into account such effect (considering then sin a ~ 1) and we interpret the magnetic 
field B as the component perpendicular to the stellar surface 

In the present work we are interested in exploring the effects due to the growth of the 
magnetic field. The growth of the magnetic field in young pulsar are not a new idea, but 
the impact on the braking index and the early dynamic of the pulsar have not been fully 
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explored. 


1 Growing magnetic fields in young pulsars 

A newborn neutron star may be exposed to a hyperaccretion phase few moments after the 
supernova explosion that originated it. The paradigm is as follows: when the core-collapse 
supernova event take place the shock is still pushing its way through the outer layers of the 
progenitor, and if it encounters a density discontinuity, a reverse shock may be generated. 
Depending on its strength and on how far out it was generated, this reverse shock can in¬ 
duce strong accretion onto the newborn neutron star on a timescales of hours. Hypercritical 
accretion results (rh > rriEdd, where friEdd is the Eddington accretion rate), in which the 
photons are trapped within the accretion flow and the energy liberated by the accretion is 
lost through neutrino emission close to the neutron star surface. After the reverse shock hits 
the neutron star surface and rebounds, a third shock develops and starts moving outward 
against the infalling matter. Once this accretion shock stabilizes it will separate the infalling 
matter from an extended envelope in quasi-hydrostatic equilibrium. Ifl6ll argued in favor 
of such scenario of late accretion onto newborn neutron stars inside supemovae and devel¬ 
oped an analytical model for the hypercritical regime. In such model, the neutrino cooling 
plays an important role in the formation of a quasi-hydrostatic envelope around the com¬ 
pact remnant. IfTTl highlighted the physical conditions presents in the formation of neutron 
stars inside supernovae: convective envelope, hyperaccretion of material and submergence 
of the magnetic field on the stellar crust. With these suggestions, llT8l presented simple ID 
ideal MHD simulations of the effect of this post-supernova hypercritical accretion on the 
newborn neutron star to show that such magnetic field submergence could occur. The result 
was a rapid burial of the magnetic field into the neutron star crust. Notice that an accreted 
mass of ~ 0.001 M 0 is enough to submerge the magnetic field (which corresponds to an 
accretion rate of ~ 100 M 0 yr -1 in few hours). It is widely accepted that the origin of 
neutron star magnetic fields is still an unsolved problem f [fT9ll . Il20ll ). Until now, two main 
mechanisms are still competing: (a) a fossil field from the progenitor compressed during 
the core collapse (magnetic flux conservation and diamagnetism of MHD turbulence), and 
(b) a proto-neutron star dynamo (generated by precollapse or the short-lived postcollapse). 
Both models are used to explain the large variety of observed field strengths. However, cur¬ 
rently is accepting the idea of magnetic field submergence by hyperacretion onto newborn 
neutron stars and subsequent reemergence, because it explains better why a neutron star at 
birth has a very low surface magnetization. In these scenarios, the magnetic field genera¬ 
tion and/or adjustment process terminates few seconds after the neutron stars birth ( IfTTlI ). 
Recently, flU; ifTOl : performing 2D-3D MHD simulations with more refined detailed physi¬ 
cal ingredients, showed that the magnetic field is submerged in the stellar surface regardless 
of their initial configuration or its strength. Although many observed pulsars show clear 
evidence of strong magnetic fields (from radio pulsar with 10 12 G to magnetars with 10 15 
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G), lower magnetic fields are, however, found in millisecond pulsars f lf2TTl ) and in neutron 
stars in low-mass X-ray binaries ( ll22l ). In such cases, hyperaccretion is thought to be the 
cause of the magnetic field reduction. In addition, there is a small group of neutron stars, 
found in young supernova remnants which exhibit little or no evidence for the presence of 
a magnetic field ( lfl8i : lf23l ). 

After hyperaccretion stopped, the magnetic field could diffuse back to the surface and 
result in a delayed switch-on of a pulsar ( lf24ll : ifTTlO . Depending on the amount of accreted 
matter, the submergence could be so deep that the neutron star may appear and remain un¬ 
magnetized for several centuries or millennia ( ||T8Tl ). This scenario was recently revisited 
by l[23ll and lfTTl and applied to study the field evolution of the CCOs (Central Compact 
Object; Il25l0 defined as X-ray sources with thermal-like spectra observed close to the cen¬ 
ters of SNRs without any counterparts in radio and gamma wavebands. With blackbody 
temperatures of about a few hundred eV and luminosities in the range 10 33 — 10 34 erg s _1 , 
they present no evidence of a pulsar magnetic field. Several of such sources are known, 
including RCW103, Cas A, Pup A, and Kes 79 ( ll26ln . In the case of magnetic fiel submer¬ 
gence, the numerical simulations were done using an ideal MHD scheme. This is possible 
because of the short duration of initial transient and due to the violence of the hypercriti¬ 
cal regime. In contradistinction, the back-diffusion of the field after hyperaccretion stops 
is due to the finite electrical conductivity of the neutron star crust matter. As a result, 
the back-diffusion becomes dependent on the thermal evolution of the star because of the 
temperature dependence of the electrical conductivity. flTl and [fTT| found that when the 
hypercritical phase is over, the magnetic field can re-emerge by a diffusion process. The 
difusion time in such case is tb ~ 10 2 3 yr. This is the timescale for grow a magnetic field 
from a low value B ~ 10 8-9 G to a high value B ~ 10 12-13 G, which depend of the early 
thermomagnetic history of the pulsar, the amount of accreted matter, the initial magnetic 
flux distribution and the electric resistivity in the crust. On the other hand, ll27ll admit that 
pulsars are rotating magnetic dipoles with growth of the magnetic field, but they did not 
propose any intrinsic mechanism for such growth. With such assumption and based on 
data from the seven pulsars with known braking indices they predict ranges for the braking 
indices of other pulsars. 

Following these ideas, it is possible to study analytically, in a first approach, the growth 
of the magnetic field when it re-emerges from the neutron star crust, and to follow its 
consequences on the pulsar spin-down. 

1.1 The modified canonical model 

The magnetic field re-diffusion from the new crust (formed in the hypercritical regime) to 
the stellar surface, under a variety of complex physical conditions, is a problem that gen¬ 
erally requires numerical modeling ( llT8ll : IfTTlO . However, valuable insight can be gained 
from purely analytical solutions, which exhibit the main properties of the spin-down dy¬ 
namics and may be applicable to various scenarios where the pulsar magnetic field is weak. 
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We think that this study is academically important because it allows us to explore the be¬ 
havior of certain parameters of the pulsar, in the early phase, under a different approach 
and without need for intensive numerical work. 

If the magnetic field has a temporal dependence through a growth function /(f), then it 
is possible to express the magnetic field as B(t ) = Bf(t), where B is the maximum mag¬ 
netic field achieved in the saturation regime. This temporal dependence affect the function 
of torque k, which is considered constant in the simple dipolar oblique rotator model. A 
similar approach was proposed initially by [f28il who proposed to replace k by some /(f) 
and them attempt to constrain /(f) using the observed value of n, but they proposed no 
growth function because his approach was focused in the thermomagnetic instabilities in 
the neutron star crust. In the present case we are interested in to study the magnetic field 
growth due a ohmic diffusion process in the stellar crust when the hypercritical phase is 
over. In such case, the structural factor of the pulsar is time-depending through the mag¬ 
netic field as, k(B(t )) = k(t), modifying the general law for the pulsar spin-down, 


Ci = -k(t)n n , k{t) = kf(t ) (4) 

Cl = —k(t)Q n , k(t ) = kf(t) (5) 

where /(f) is an analytical function that allows the growth of the magnetic field from an 
low initial value to a maximum value in the saturation regime, and k is the standard constant 
of the pulsar. Changes in the moment of inertia are neglected in this approach (see next 
section in which such changes are investigated). Performing a similar treatment as in the 
canonical model, we find the relationships for the overall properties of the young pulsar 
with growth of the magnetic field: characteristic age, braking index, spin-down luminosity 
and period. 

The characteristic age is found integrating the modified general law <[5]) as, 


1 

m 


t 

T o + J f(t)dt 

o 


( 6 ) 


Notice that, depending on the function /(f) proposed, the characteristic age changes sub¬ 
stantially from its canonical value. 

The braking index n is found by a straightforward differentiation of the modified gen¬ 
eral law ([5]) as, 


/(f) If /(f) P 

n = n* + tttt — = n* - 


f(t) Cl 


f(t)P ’ 


(7) 


where the theoretical braking index is n* = 3 and /(f) = (If /(It is the temporal derivative 
of /(f). Notice that to obtain n < 3 is required magnetic field growth through /(f) > 0, 
because Cl < 0 always (the pulsar is slowing). 
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In addition, in the modified canonical model, the spin-down luminosity and the period 
are given by, 
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In this case, the temporal evolution of the spin-down luminosity and rotation period look 
more complex that in the canonical model, and depend exclusively of the growth function 

fit)- 


1.2 The growth function fit) 

Regardless of which analytical function /(f) is proposed, one of its main requirements is 
that for t > t b , the canonical properties for the dipolar oblique rotator model (character¬ 
istic age, braking index, spin-down luminosity and period) are recovered. That is, when 
the magnetic field reaches the saturation regime, the evolution of the pulsar follows the 
canonical model. Thus, it is required that /(f) allows the growth of the magnetic field 
rapidly, in a diffusion timescale t b , from a low value to a canonical value. As already 
mentioned, magnetic field strengths inferred from the observed pulsar population range be¬ 
tween B ~ 10 8 G for recycled (or millisecond) pulsar, B ~ 10 12 G for normal pulsars up to 
B ~ 10 15 G for magnetars. Therefore, the growth function must satisfy that /(f — > 0) = e 
and /(f > t b ) = 1, where e = B 0 /B lisa parameter related with the magnetic field 
strength B, in the saturation regime, and B 0 is the initial low magnetic field. An expo¬ 
nential growth function is most appropriate in this case, which satisfies the aforementioned 
requirements, 


/(f) — e + 


1 — exp 



(9) 


The parameters e and t b involve the most relevant physical ingredients of the problem: 
magnetic fields, resistivity in the crust, MHD instabilities, submergence depths and more. 
However, the oversimplified growth function ([9]) match very well with those numerical 
magnetic curve obtained by [|T8l & ffTTil for standard parameters. Fig.[2]shows the evolution 
of the magnetic field through the growth function /(f) and the characteristic age evolution 
for young pulsars. In (a) we show B(t) for several values of e and for a fixed t b = 10 3 yr, to 
analyze its effect on the growth of the magnetic field. Notice the rapid magnetization of the 
pulsar. Due to the constriction on t b , the magnetic field is saturated more rapidly for low 
values of e and slower for higher values of the same parameter. That is, the magnetic field 
is almost constant during a time and may take longer to start the growth process for high 
values of e, while for low values of e it may take less time to start such process of growth 
(from B = 10 5 G to B = 10 12 G). But in any case, the growth of the magnetic field is very 
efficient (occurring in a timescale t b ). In (b) we show the characteristic age r for several 
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(a) B(t ) for several e, fixing r B 



Figure 2: (color online) In (a) we show the growth of magnetic field through function /(£), 
in logarithmic scale, varying e for a fixed r B = 10 3 yr. We show B(t ) for several values 
of e. The canonical model, B = 10 12 G (or equivalently f(t ) = 1), is superimposed for 
comparison purposes. In (b) we show the modified characteristic age of the pulsar for a 
fixed value of e = 0.001 and varying the diffusion timescale r B . The canonical parameter 
t is superimposed for comparison purposes. 
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(a) E for several tb, fixing e 



(b) n vs age for several tb, fixing e 

Figure 3: (color online) In (a) we show the modified spin-down luminosity of the pulsar 
for a fixed value of e = 0.001 and varying the diffusion timescale tb- In (b) we show the 
modified braking index of the pulsar for same parameters as (a). The canonical parameters, 
E and n, are superimposed for comparison purposes. 
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values of tb and for a fixed e = 0.001. Notice that initially the modified r is different from 
the canonical value, but after a diffusion timescale tb it evolve like the familiar canonical 
model. The main reason for this behavior is the parameter e in the growth function f(t ) 
because when t —> 0, r —* r 0 /e, being r 0 the canonical initial spin-down timescale of the 
pulsar. That is, for an initial low magnetic field B 0 = 10 9 G, the corresponding initial 
characteristic age is three orders of magnitude higher that the canonical case, which is 
chosen as r 0 = 500 yr for this example, in particular. 

In Fig. [3] we show both the spin-down luminosity and the braking index evolution. In 
(a) the modified spin-down luminosity is shown for several values of tb and for a fixed 
e = 0.001. Notice that initially the modified spin-down luminosity E is different from 
the canonical value, but after a diffusion timescale tb it evolve like the familiar canonical 
model. Again, the reason for this behavior is the parameter e in the growth function /(f) 
because when t —> 0, E —>■ eE 0 , being E 0 the canonical initial spin-down luminosity of 
the pulsar. That is, for an initial low magnetic field B 0 = 10 9 G the corresponding initial 
spin-down luminosity is three orders of magnitude lower that the canonical case, which is 
chosen as E {) = 10 38 erg s -1 for the present example. In (b) we show the modified braking 
index n for several values of tb and for a fixed e = 0.001. Notice that n depend of tb and 
e, but has a direct dependence with the period P and its derivative P too. In this case, when 
t —> 0, n —> 3 — PP~ 1 /eTf-j. The plot shows the behavior of the braking index for standard 
values of period and derivative of period of pulsars. Again, after a diffusion timescale tb 
it evolve like the familiar canonical model. That is, the growth of the magnetic field allow 
the braking index be less than the canonical value. After a timescale t > tb, the braking 
index parameter tends to n = 3. 


2 Conclusions 

We have revisited the magnetic torque problem from two different approaches: changes 
due growth of the magnetic field in the pulsar and changes in the moment of inertia. In 
the first approach, we studied the early evolution of the magnetic field through a growth 
function /(f), which has a exponential behavior, and followed its consequences on the 
overall properties of RPP, focusing on very young pulsars. The main motivation for the 
present study is due to the recent observations of a population of young pulsars, showing 
no magnetic field which is characteristic of these objects ( [12311 . ll27lD . It has already been 
shown in previous work that the lack of magnetic field in young neutron star may be due to 
the fact that it is submerged into stellar crust during an early hypercritical phase. 

However, it is possible that when the hypercritical phase ends, this magnetic field may 
emerge to the surface, due to MHD instabilities or buoyancy effects, and be displayed as 
a delayed pulsar. We find that the growth of the magnetic field change the early dynamics 
of the pulsar and modifies the RPP canonical model in the early regime. In particular, we 
observe that initially the characteristic age is higher than expected, depending on how the 
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magnetic field grows. This implies that the pulsar looks older than it actually is. As was 
suggested by ifTSTl . one way to identify such a pulsars would be through an association with 
a supernova remnant, which is a challenging task. Also, the initial spin-down luminosity is 
lower than that expected from the canonical model, however, it eventually increases until 
reaching the canonical behavior (again due to the growth of the magnetic field). This is 
interesting because if the surface magnetic field is low (in the early history of the pulsar), 
the system radiates its energy inefficiently, allowing a low initial spin-down luminosity 
and a high initial characteristic age. In contrast, the spin period evolves in similar way 
to the canonical model, due to its weak dependence on /(f). Thus it remains practically 
unchanged for t < tb and a slight change is observed during the growth of B. 

Finally, we show that the braking index is less than the expected value from the theo¬ 
retical canonical model, n = 3, when a growth of magnetic field process is present. We 
believe that such behavior may help explain (even if only partially) the low values observed 
in the few young pulsar where this parameter has been measured accurately. Nevertheless, 
very accurate values of the age of the supernova remnant are necessary for a better esti¬ 
mate of diffusion timescale tb- Of course, we are assuming that the RPP model is correct 
in general, and we augment it by including the possibility of magnetic field growth. For 
instance, the Crab pulsar has, until now, the better initial spin determination, which was ob¬ 
tained independently: P 0 = 19 ms. The values of period and its derivative (P = 33 ms and 
P_i 4 = 42 s s -1 ) obtained from timing observation allow to estimate the current braking 
index as n = 2.51 which tends very slowly to its canonical value. The current age is known 
from historical records (SN1054) and it is f ~ 960 yr. With these parameters it is possible 
to estimate the rotational energy as, E rot ~ 2 x 10 46 / 45 P _ 2 (s)erg ~ 2 x 10 49 erg. The spin- 
down luminosity is given by, E rot ~ —4 x 10 32 / 45 P” 3 (s)P_i 4 erg s " 1 ~ — 5 x 10 38 erg s _1 . 
The inferred superficial magnetic field is then, B 12 — 6 (/ 45 P (s) P_ 1 4 ) 0 ’ 5 / P;' sin a G ~ 7, 
where Z 45 is the moment of inertia in units of 10 45 g cm 2 and P _ i4 is the period derivative 
in units of 10 ~ 14 s s' 1 . Also, if the magnetic field grew three orders of magnitude in such 
system (e = 0 . 001 ), the estimated diffusion timescale is given by, t b = 520 yr. ifTTl and 
IfUll gave estimates of the difussion timescale performing 1D/2D numerical simulations. 
Meanwhile, we do a phenomenological study of the behavior of the magnetic field and 
studies other properties of the pulsar. We propose that it would be possible to estimate this 
parameter more accurately using the observed data for each pulsar with such simulations. 
This will accurately calculate the expected values of n showing that the growth of the mag¬ 
netic field in these pulsars is feasible. This approach is prudent and necessary to carry out 
but is beyond the scope of this work. 
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